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Abstract 

The compelling experimental evidences for oscillations of solar and atmospheric 
neutrinos imply the existence of 3-neutrino mixing in vacuum. We briefly review the 
phenomenology of 3-is mixing, and the current data on the 3-neutrino mixing parame- 
ters. The open questions and the main goals of future research in the field of neutrino 
mixing and oscillations are outlined. The predictions for the effective Majorana mass 
|<m>| in {(3(3)q v — decay in the case of 3-v mixing and massive Majorana neutrinos are 
reviewed. The physics potential of the experiments, searching for (/3(3)o u — decay and 
having sensitivity to |<m>| ^ 0.01 eV, for providing information on the type of the 
neutrino mass spectrum, on the absolute scale of neutrino masses and on the Majorana 
CP-violation phases in the PMNS neutrino mixing matrix, is discussed. 
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■ 1 Introduction 

There has been a remarkable progress in the studies of neutrino oscillations in the last 
several years. The experiments with solar, atmospheric and reactor neutrinos PH 121 El IH EH 
Ej have provided compelling evidences for the existence of neutrino oscillations driven by 
\q . nonzero neutrino masses and neutrino mixing. Evidences for oscillations of neutrinos were 
obtained also in the first long baseline accelerator neutrino experiment K2K |Hj. 

The idea of neutrino oscillations was formulated in It was predicted in 1967 [10J 
that the existence of v e oscillations would cause a "disappearance" of solar v e on the way 
to the Earth. The hypothesis of solar v e oscillations, which (in one variety or another) were 
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considered from ~1970 on as the most natural explanation of the observed PP solar u e deficit 
(see, e.g., refs. fH\ IT2"1 IT4"] ) . has received a convincing confirmation from the measurement 
of the solar neutrino flux through the neutral current reaction on deuterium by the SNO 
experiment [Sj, and by the first results of the KamLAND (KL) experiment [Zj. The combined 
analysis of the solar neutrino data obtained by Homestake, SAGE, GALLEX/GNO, Super- 
^ . Kamiokande (SK) and SNO experiments, and of the KL reactor v e data |7j, established the 
large mixing angle (LMA) MSW oscillations/transitions ^2] as the dominant mechanism 
at the origin of the observed solar v e deficit (see, e.g., JH]). The Kamiokande experiment 
16J provided the first evidences for oscillations of atmospheric and u^, while the data 
of the SK experiment made the case of atmospheric neutrino oscillations convincing [HI E] . 
Indications for i/-oscillations were reported by the LSND collaboration [T7j . 

The latest contributions to these magnificent progress are the new SK data on the L/E- 
dependence of the /i-like atmospheric neutrino events [18 , L and E being the distance 
traveled by neutrinos and the neutrino energy, and the new spectrum data of the KL and 
K2K experiments ^01 120] • For the first time the data exhibit directly the effects of the 

*Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 Sofia, 
Bulgaria. 
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oscillatory dependence on L/E and E of the probabilities of z/-oscillations in vacuum [2JJ- As 
a result of these developments, the oscillations of solar u e , atmospheric and accelerator 
(at L ~250 km) and reactor v e (at L ~180 km), driven by nonzero iz-masses and z/-mixing, 
can be considered as practically established. 

2 The Neutrino Mixing Parameters and ((3(3)o u — Decay 

The SK atmospheric neutrino and K2K data are best described in terms of dominant 2- 
neutrino — > u T {y^ — > u T ) vacuum oscillations. The best fit values and the 99.73% C.L. 
allowed ranges of the atmospheric neutrino oscillation parameters read 

I Am| I = 2.1 x 1(T 3 eV 2 , sin 2 29 A = 1.0 ; 

(1) 

| Ami I = (1-3 - 4.2) x 1(T 3 eV 2 , sin 2 29 A > 0.85. 

It should be noted that the signs of Am\ and of cos26'a, if sin 2 29a ^ 1.0, cannot be 
determined using the existing data. 

Combined 2-neutrino oscillation analyses of the solar neutrino and the new KL spectrum 
data show ^3 1221 that the ^©-oscillation parameters lie in the low-LMA region : Am^ = 
(7.9±gJ) x 10" 5 eV 2 , tan 2 9 & = (0.40ig;g?). The high-LMA solution (see, e.g., [HJ) is excluded 
at ~ 3.3cr. Maximal z/ Q -mixing is ruled out at ~ 6a; at 95% C.L. one finds cos 29 Q > 0.28 [2~2*] . 
which has important implications (see further). One also has: Am| /|Am A | ~ 0.04 <C 1. 

The evidences for i/-oscillations obtained in the solar and atmospheric neutrino and KL 
and K2K experiments imply the existence of 3-v mixing in the weak charged lepton current: 
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f& = Uij u jL , I = e, /i, r, (2) 



3=1 



where vn are the flavour neutrino fields, i/jl is the field of neutrino Uj having a mass rrij 
and U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix 0123, U = t/pMNS- 
All existing i/-oscillation data, except the data of LSND experiment J7j, can be described 
assuming 3-v mixing in vacuum and we will consider this possibility in what follows 1 . 

The PMNS matrix can be parametrized by 3 angles, and, depending on whether the 
massive neutrinos Uj are Dirac or Majorana particles, by 1 or 3 CP-violation (CPV) phases 
[2Ell2Ij. In the standardly used parameterization (see, e.g., |2Hj), C^pmns has the form: 



PMNS 



/ C12C13 S12C13 S13 \ 

-S12C23 - Ci 2 S23Sl 3 e i ' 5 C12C23 - Si2S23Sl3e i<5 s 2 3Ci 3 e J<5 
V S12S23 — Ci2C23Si3e l<5 — C12S23 — Si2C23Si3e l<5 Cl^Vi^ ) 



diag(l, e l 2 1 ; e 1 ^) 



(3) 

where c^- = cos 6^-, Sy = sin^, the angles 9^ = [0, 7r/2], 5 = [0, 2tt] is the Dirac CPV phase 
and «2i, «3i are two Majorana CPV phases |26| 12 7j. One can identify Am^ = Am| x > 0. 
In this case |Aml | = (Am^l = |Am 2 2 | Am^, #12 = 9 Q , 6*23 = #a- The angle #13 is lim- 
ited by the data from the CHOOZ and Palo Verde experiments [2H]. The presently existing 

1 In the LSND experiment indications for — > v e oscillations with (Am 2 )LSND — 1 eV 2 were obtained. 
The minimal 4-v mixing scheme which could incorporate the LSND indications for oscillations is strongly 
disfavored by the data [21] ■ The ^-oscillation explanation of the LSND results is possible assuming 5-v 
mixing The LSND results are being tested in the MiniBooNE experiment |25| . 
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atmospheric neutrino data is essentially insensitive to 9\s satisfying the CHOOZ limit jSj- 
The probabilities of survival of solar u e and reactor i> e , relevant for the interpretation of the 



solar neutrino, KL and CHOOZ data, depend in the case of interest, 
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where P 2 " is the solar v e survival probability jHUj corresponding to 2-v oscillations driven 
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in which the solar e number density iV e is replaced by N e cos # 
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P 2u + Pq" osc , Pq osc being an oscillating term [3U| and 



p2v 



\-P') cos 29™ (t ) cos 20i2, 



-27T7-0 n tT 1 



(4) 
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Here P^ is the average probability |32[l30j. P' is the "double exponential" jump probability 
30J and r$ is the "running" scale-height of the change of N e along the //-trajectory in the 
Sun 2 mj E3 El- In the LMA solution region P^ osc 9* [H]. Using the expressions for 
Pkl? Pchooz an d -f© 1 ' gi yen above, the 3a allowed range of \Am\ | from [3], and performing 
a combined analysis of the solar neutrino, CHOOZ and KL data, one finds [22] : 



sin 



9 13 < 0.055, 99.73% C.L. 



(5) 



Similar constraint is obtained from a global 3-z/ oscillation analysis of solar, atmospheric and 
reactor neutrino data jSHEEj- A combined 3-z/ oscillation analysis of the solar neutrino, 
CHOOZ and KL data shows also [22J that for sin 2 #13 ^ 0.02 the allowed ranges of Am\\ 
and sin 2 6* 2 i do not differ substantially from those derived in the 2-v oscillation analyzes (see, 
e.g., ref. [H]). The best fit values, e.g., read 3 



Am 2 ! 



e.g. 

8.0 x 10~ 5 eV 2 



sin 2 9 
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0.28. 



(6) 

In Fig. ^we show the allowed regions in the Am\\ -sin 2 9\ 2 plane, obtained in a 3-z/ oscillation 
analysis of the solar neutrino, KL and CHOOZ data for few fixed values of sin 2 # 13 . 

As we have seen, the fundamental parameters characterizing the 3-neutrino mixing are: 
i) the 3 angles #12, #23, #13, h) depending on the nature of Uj - 1 Dirac (5), or 1 Dirac + 
2 Majorana (£,0:21,031), CPV phases, and iii) the 3 neutrino masses, mi, m 2 , m 3 . It is 
convenient to express the two larger masses in terms of the third mass and the measured 
Am 2 & =Am\ x >0 and Am\ . We have remarked earlier that the atmospheric neutrino 



and K2K data do not allow one to determine the sign of Am A . This implies that if we 
identify Am\ with Am^^ in the case of 3-neutrino mixing, one can have Am^j > or 

2 The analyses and the extensive numerical studies performed in |.'S."SI LI lj show that expression (0J for P^ 
provides a high precision description of the average solar v e survival probability in the Sun for any values 
of Am 2 ^ and 6i2 (the relevant error does not exceed ~(2-3)%), including the values from the LMA region. 
The results obtained recently in (HSJ imply actually that the use of the double exponential expression for P' 
for description of the LMA transitions brings an imprecision in Pq which does not exceed ~ 10~ 6 . 

3 The best fit value of sin 2 #13 = 0.004 is different from zero but not at statistically significant level 
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Figure 1: The 90%, 95%, 99% and 99.73% C.L. allowed regions in the Am 2 21 -sin 2 12 plane, 
obtained in a 3-v oscillation analysis of the solar neutrino, KL and CHOOZ data |2"2"] . 



^ m 3i(2) < 0. The two possible signs of Am\ correspond to two types of z/-mass spectrum: 



with normal hierarchy, mi < m 2 < m 3 , Am\ =Am| 1 > 0, m 2 ( 3 ) = {m\ + Am 21 
with inverted hierarchy 4 , m 3 < mi < m 2 , Am\ =Am§ 2 <0, m 2 =(m| — Am| 2 )^, etc. 



(31) ;2, and 
i 



i 

I 2 rvj 



0.009 eV, m^\Am\ 



The neutrino mass spectrum can also be 

- Normal Hierarchical (NH): mi< m 2 ^Cm 3 , m 2 =(AmQ 

- Inverted Hierarchical (IH): m 3 m\ < m 2 , with mi j2 = |Am| | s ~0.045 eV; or 

- Quasi- Degenerate (QD): mi = m 2 = m 3 = m , m 2 3> |Am^ |, m ^ 0.20 eV. 

Neutrino oscillation experiments allow to determine differences of squares of neutrino 
masses, but not the absolute values of the masses, or rninirrij). Information on the absolute 
scale of v- masses can be derived in 3 H /3-decay experiments [SHI EDI E] and from cosmological 
and astrophysical data (see, e.g., ref. |12!)- The currently existing most stringent upper 

and Mainz [3T] experiments: 



bounds on the v e mass were obtained in the Troitzk 



m De < 2.3 eV (95% C.L). 



(7) 



We have m Pe = mi j2j3 in the case of QD z/-mass spectrum. The KATRIN experiment ^T] is 
planned to reach a sensitivity to m Pe ~ 0.20 eV, i.e., to probe the region of the QD spectrum. 
The CMB data of the WMAP experiment were used to obtain the upper limit [4*3*j : 

2 m j < (0.7-2.0) eV (95% C.L.), (8) 

j 

where we have included a conservative estimate of the uncertainty in the upper limit (see, 
e.g., |42j). The WMAP and future PLANCK experiments can be sensitive to J2j m j — 0.4 



4 In the convention we use (called A), the neutrino masses are not ordered in magnitude according to their 
index number: Am^ < corresponds to TO3 < mi < 7712. We can also always number the neutrinos with 
definite mass in such a way that [37| m\ < mi < 777,3. in this convention (called B), we have in the case of 
inverted hierarchy spectrum: AttZq =A777§ 2 , Am^ =Am§ 1 . Convention B is used, e.g., in [2HlEHj- 
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eV. Data on weak lensing of galaxies by large scale structure, combined with data from the 
WMAP and PLANCK experiments may allow one to determine (m\ + m 2 + m 3 ) with an 
uncertainty of 5 ~ (0.04 — 0.10) eV (see [12] and the references quoted therein). 

The type of neutrino mass spectrum, i.e., sjnfAmjj), can be determined by studying 
oscillations of neutrinos and antineutrinos, say, <-+ v e and <-> // e , in which matter 
effects are sufficiently large. This can be done in long base-line //-oscillation experiments 
(see, e.g., |Bj). If sin 2 2#i 3 ^0.05 and sin 2 6*23 ^0.50, information on sgn(/S.m\ 1 ) might be 
obtained in atmospheric neutrino experiments by investigating the effects of the subdominant 
transitions v^ e ) u e{^) and v^ e ) —> Ve(p) °f atmospheric neutrinos which traverse the Earth 
For Pfj,< e ) {or // M ( e )) crossing the Earth core, new type of resonance-like enhancement of 
the indicated transitions takes place due to the (Earth) mantle- core constructive interference 
effect (neutrino oscillation length resonance (NOLR)) 42] 5 . As a consequence of this effect 6 
the corresponding v^u) (orv^) transition probabilities can be maximal |47j . For Am^ > 0, 
the neutrino transitions zv e ) — * v e {^) are enhanced, while for Am^ < the enhancement of 
antineutrino transitions — > u e ^ takes place, which might allow to determine s^n(Am| 1 ). 

After the spectacular experimental progress made in the studies of //-oscillations, further 
understanding of the structure of neutrino mixing, of its origins and of the status of CP- 
symmetry in the lepton sector, requires a large and challenging program of research to be 
pursued in neutrino physics. The main goals of this research program should include [30] : 

- High precision measurement of neutrino mixing parameters which control the solar and 

the dominant atmospheric neutrino oscillations, Am^, #21, and Am^, 6*23- 

- Measurement of, or improving by at least a factor of (5 - 10) the existing upper limit on, 

813 - the only small mixing angle in £/pmns- 

- Determination of the sgn(Am\^) and of the type of //-mass spectrum (NH, IH, QD, etc.). 

- Determination or obtaining significant constraints on the absolute scale of //-masses. 

- Determination of the nature-Dirac or Majorana, of massive neutrinos Uj. 

- Establishing whether the CP-symmetry is violated in the lepton sector a) due to the Dirac 

phase 5, and/or b) due to the Majorana phases a^i and a^i if Vj are Majorana particles. 

- Searching with increased sensitivity for possible manifestations, other than flavour neutrino 

oscillations, of the non-conservation of the individual lepton charges Li, I = e, //, r, such 
as [A — > e + 7, r — * fi + 7, etc. decays. 

- Understanding at fundamental level the mechanism giving rise to neutrino masses and 
mixing and to L;— non- conservation. This includes understanding the origin of the patterns 
of //-mixing and //-masses suggested by the data. Are the observed patterns of //-mixing and 
of Am^ 31 related to the existence of new symmetry of particle interactions? Is there any 
relations between quark mixing and neutrino mixing? Is #23 = tt/4, or #23 > 7r/4 or else 
#23 < vr/4? Is there any correlation between the values of CPV phases and of mixing angles 
in t/pMNS? Progress in the theory of neutrino mixing might also lead, in particular, to a 
better understanding of the origin of baryon asymmetry of the Universe |57j . 

5 For the precise conditions of the mantle-core (NOLR) enhancement see |4fil 147) . 

6 The Earth mantle-core (NOLR) enhancement of neutrino transitions differs @B] from the MSW one. 
It also differs |4(j| 147] from the parametric resonance mechanisms of enhancement discussed in |48j : the 
conditions of enhancement found in 0S| are not realized for the neutrino transitions in the Earth. In it 
is erroneously concluded that the ^V e ) v e{n) transitions of atmospheric neutrinos crossing the Earth core 
cannot be enhanced by the interplay of the transitions in the Earth mantle and of those in the Earth core. 
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The mixing angles, 8 2 i, #23 and #13, the Dirac CPV phase 5 and Am^ and Arn^j can, 
in principle, be measured with a sufficiently high precision in //-oscillation experiments (see, 
e.g., |H| I51j). These experiments, however, cannot provide information on the z/-mass scale 
and on the nature of massive neutrinos Vj, they are insensitive to the Majorana CPV phases 
«2i,3i 15 2 j . Establishing whether Uj have distinct antiparticles (Dirac fermions) or not 
(Majorana fermions) is of fundamental importance for understanding the underlying sym- 
metries of particle interactions and the origin of //-masses. If Vj are Majorana fermions, 
getting experimental information about the Majorana CPVphases in C/pmns would be a re- 
markably challenging problem [HZ1 EH IHH1 EH ISHj • The phases 0121,31 can affect significantly 
the predictions for the rates of (LFV) decays \i — >• e + 7, r— > + 7, etc. in a large class 
of supersymmetric theories with see-saw mechanism of //-mass generation (see, e.g., 56 ). 
Majorana CPV phases might be at the origin of baryon asymmetry of the Universe |57j - 

In the present article we will review the potential contribution the studies of neutrinoless 
double beta ((f3f3)o u — ) decay of even-even nuclei, (A, Z) — ► (A, Z + 2) + e~ + e~, can make 
to the program of research outlined above. The (/3/3)o!,-decay is allowed if the neutrinos with 
definite mass are Majorana particles (for reviews see, e.g., [131 EHl EHl ED])- Let us recall that 
the nature - Dirac or Majorana, of the massive neutrinos Uj, is related to the symmetries of 
particle interactions. Neutrinos Vj will be Dirac fermions if the particle interactions conserve 
some lepton charge, e.g., the total lepton charge L. If there does not exist any conserved 
lepton charge, neutrinos with definite mass can be Majorana particles. As is well-known, 
the massive neutrinos are predicted to be of Majorana nature by the see-saw mechanism of 
neutrino mass generation (see, e.g., [HI]), which also provides a very attractive explanation 
of the smallness of the neutrino masses and - through the leptogenesis theory [57j, of the 
observed baryon asymmetry of the Universe. 

If the massive neutrinos Uj are Majorana fermions, processes in which the total lepton 
charge L is not conserved and changes by two units, such as K + — ► 7r _ + fi + + fi + , fi + + 
(A, Z) — > (A, Z + 2) + etc., should exist. The process most sensitive to the possible 
Majorana nature of neutrinos Uj is the (PP)q u — decay (see, e.g., O If the (Pf3)ov- 

decay is generated only by the (V-A) charged current weak interaction via the exchange 
of the three Majorana neutrinos Vj (rrij ^ 1 eV), which will be assumed in what follows, 
the dependence of the (/3/3)o^-decay amplitude A((3/3)o u on the neutrino mass and mixing 
parameters factorizes in the effective Majorana mass <m> (see, e.g., l59|): 

A(PP)ov ~ <m> M , (9) 
where M is the corresponding nuclear matrix element (NME) and 

\<m>\ = \ mi \U el \ 2 + m 2 \U e2 \ 2 e ia21 +m 3 \U e3 \ 2 e ia31 \ . (10) 

If CP-invariance holds 7 , one has 63 a 2 \ = kit, a 31 = k'n, where k, k' = 0, 1, 2, and 

77 21 = e*° 21 = ±1, 7731 = e"* 31 = ±1 (11) 

represent the relative CP-parities of Majorana neutrinos V\ and u 2 , and V\ and z/ 3 , respec- 
tively. It follows from eq. (jlOj) that the measurement of |<m>| will provide information, in 
particular, on rrij. As eq. indicates, the observation of (/3/3)ojy-decay of a given nucleus 

7 We assume that rrij > and that the fields of the Majorana neutrinos Vj satisfy the Majorana condition: 



j 

C(i?j) T = vj, j — 1, 2, 3, where C is the charge conjugation matrix 



6 



and the measurement of the corresponding half-life, would allow to determine |<m>| only 
if the value of the relevant NME is known with a relatively small uncertainty. 

The experimental searches for (fl(3)ov— decay have a long history (see, e.g., [2E1EH])- The 
best sensitivity was achieved in the Heidelberg- Moscow 76 Ge experiment |64j : 

|<m>| < (0.35 - 1.05) eV, 90% C.L. (12) 

where a factor of 3 uncertainty associated with the calculation of the relevant nuclear matrix 
element is taken into account. A positive signal at >3cr, corresponding to |<m>| = 
(0.1 — 0.9) eV at 99.73% C.L., is claimed to be observed in (HSj- This result will be checked 
in the currently running and future (/3/3)o^-decay experiments. However, it may take a long 
time before comprehensive checks could be completed. Two experiments, NEM03 (with 
100 Mo and 82 Se) jHi and CUORICINO (with 130 Te) jHT], designed to reach sensitivity to 
|<m>| ~(0.2-0.3) eV, are taking data. Their first results read, respectively (90% C.L.): 

|<m>| < (0.7-1.2) eV |5E|, |<m>| < (0.3 - 1.6) eV (13) 
where estimated uncertainties in the NME are accounted for. A number of projects aim to 
reach sensitivity to |<m>| -(0.01-0.05) eV |52E2|: CUORE ( 130 Te), GENIUS ( 76 Ge), EXO 
( 136 Xe), MAJORANA ( 76 Ge), MOON ( 100 Mo), XMASS ( 136 Xe), etc. These experiments can 
probe the region of IH and QD spectra and test the positive result claimed in [65J. 

The (/3/3) ^-decay experiments are presently the only feasible experiments capable of es- 
tablishing the Majorana nature of massive neutrinos ^SUSHIIES]. As we will discuss in what 
follows, a measurement of a nonzero value of |<m>| ^ 10~ 2 eV: 

- Can give also information on neutrino mass spectrum [HB1EB1IHE] (see also jZOj)- 

- Can provide unique information on the absolute scale of neutrino masses (see, e.g., jHHl)- 

- With additional information from other sources ( 3 H /?-decay experiments and/or cosmo- 
logical/astrophysical data) on the absolute z/-mass scale, can provide unique information on 
the Majorana CPV phases a 2 i and/or a 3 i |2H1 E3 E3 EHl E] • 

3 Properties of Majorana Neutrinos: Brief Summary 

The properties of Majorana fields (particles) are very different from those of Dirac fields 
(particles). A massive Majorana neutrino Xk can be described (in local quantum field theory) 
by 4-component complex spin 1/2 field Xk(%) which satisfies the Majorana condition: 

C (Xk(x)) T = SkXk(x), |&| 2 =1. (14) 

where C is the charge conjugation matrix. The Majorana condition is invariant under proper 
Lorentz transformations. It reduces by 2 the number of independent components in Xk(%)- 
The condition (fT4*j) is invariant with respect to U(l) global gauge transformations of the 
field Xk{%) carrying a U(l) charge Q, Xfc( x ) ~~ * e%a ®Xk{ x )i only if Q — 0. As a result, i) 
Xk cannot carry nonzero additive quantum numbers (lepton charge, etc.), and ii) the field 
Xki x ) cannot "absorb" phases 8 . Thus, Xk( x ) describes 2 spin states of a spin 1/2, absolutely 
neutral particle, which is identical with its antiparticle, Xk = Xk- If CP- invariance holds, 
Majorana neutrinos have definite CP-parity r)cp(xk) — =t^ : 

U C p Xk(x) U C p = r] CP (xk) 7o Xk(x'), r] CP (xk) = ±i ■ (15) 

8 This is the reason why the PMNS matrix contains two additional CPV phases in the case of massive 
Majorana neutrinos |26| . 
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It follows from the Majorana condition that the currents Xk{%)O l Xk{x) = 0, for O l = 7 a ; 
&a/3', &apl5- This means that Majorana fermions (neutrinos) cannot have nonzero U(l) 
charges and intrinsic magnetic and electric dipole moments. 

Finally, if $?(x) is a Dirac field and we define the standard propagator of \l/(x) as 

< 0|r(* a (x)^(y))|0 >= S^(x - y) , (16) 

one has 

< 0\T(* a (x)*p(y))\0 >= , < Q\T(V a (x)* (y))\O >= . (17) 
In contrast, a Majorana neutrino field Xk( x ) has, in addition to the standard propagator 

< o\T(xUx)xk{i(y))\Q >= O s ~ v) > ( 18 ) 

two non-trivial non-standard (Majorana) propagators 9 

< 0\T( Xka (x) Xk ^y))\0 >= -CS%{x - y)C 5p , 

< a\T(xUx)xk P (y))\o >= S cjs™(x - y ). 

This result implies that if Vj(x) in eq. (j2J) are massive Majorana neutrinos, (/?/?)o y -decay can 
proceed by exchange of virtual neutrinos Vj since < Q\T{vj a {x)vjp(y))\Q >^ 0. 

4 Predictions for the Effective Majorana Mass 

The predicted value of |<m>| depends in the case of 3 — v mixing on [72] (see also |2"H1 ITU]): 

i) Am\ = AmLgji, ii) Q = Q\i and Am| = Am^, iii) the lightest neutrino mass, 
minim f) and on iv) the mixing angle #13. In the convention (A) employed by us, one 

has \U cX \ 2 = cos 2 Ml - |t4 3 | 2 ), l^l 2 = sin 2 # (l - l^ 2 ), and |^ e3 | 2 = sin 2 13 . * 

Given Awq , Am\ , 9 Q and sin 2 ^, the value of |<m>| depends strongly on the type 
of the neutrino mass spectrum as well as on the values of the two Majorana CPV phases of 
the PMNS matrix, 0:21,31 (see eq. (fTTIjl ). In the case of QD spectrum (mi = m<i = m 3 = m , 
?71q \Am 2 A |, Ami), |<m>| is essentially independent on Am\ and Am 2 & , and the two 
possibilities, Am\ > and Am\ < 0, lead effectively to the same predictions for |<m>|. 

Normal Hierarchical Neutrino Mass Spectrum. In this case one has 



\<m>\ 



>i cos 2 9 Q + e ia21 ^Am% sin 2 6 Q ) cos 2 9 13 + \] Am\ sin 2 13 e ia31 

2 c lr, 2 fl„ n^Q 2 fl,„ -I- J Am 2 oin 2 fl, „ o*( a 31 ~"21 ) 



Am & sin 9 & cos 6\ 3 + J Am A sin #i 3 e 



(20) 
(21) 



where we have neglected the term ~ mi in eq. (|21)1. Although neutrino z/ x effectively 
"decouples" and does not contribute to |<m>|, eq. (J2TJ), the value of |<m>| still depends 
on the Majorana CPV phase CK32 = 031—021- This reflects the fact that in contrast to the case 
of massive Dirac neutrinos (or quarks), CP- violation can take place in the mixing of only two 
massive Majorana neutrinos [2E]. Further, since [221 J Am^ ^ 9.5 x 1CT 3 eV, sin 2 6 Q ^ 0.36, 



Am\ ^ 5.4 x 10 2 eV [22] (at 90% C.L.), and the largest neutrino mass enters into the 
expression for |<m>| with the factor sin 2 6 l i 3 < 0.055, the predicted value of |<m>| is 
typically ~ few x 10" 3 eV: for sin 2 13 = 0.04 (0.02) one finds |<m>| < 0.005 (0.004) 
eV. Using the best fit values of the indicated parameters (see eqs. (JTj) and ©) we get 

9 For further detailed discussion of the properties of Majorana neutrinos (fermions) see, e.g., [T3| . 
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|<m>| ^ 0.0044 (0.0035) eV. It follows from eq. (|2Uj) and the allowed ranges of values of 
AttIq , Am\ , sin 2 #0, sin 2 i3 as well as of the lightest neutrino mass m x and CPV phases 
0^21,31 that in the case of spectrum with normal hierarchy there can be a complete cancellation 
between the three terms in eq. (120 jl and one can have [HB] |<m>| = 0. 

Inverted Hierarchical Spectrum. For IH neutrino mass spectrum (see, e.g., |28j) 



m 3 <rai = m 2 = J\Am\ \ = J Am^. Neglecting 7713 sin 2 6*13 in eq. (fTUJ) . we find [37J: 



Am 2 A I cos 2 013^1 - sin 2 20© sin 2 (22) 

Even though one of the three massive Majorana neutrinos "decouples", the value of |<m>| 
depends on the Majorana CP- violating phase a 2 i- Obviously, 



Ami I I cos 20 | cos 2 #13 < |<m>| < \j\Am\ | cos 2 #13. (23) 

The upper and the lower limits correspond respectively to the CP-conserving cases. Most 
remarkably, since according to the solar neutrino and KamLAND data cos 29 & ~ 0.40, we 
get a significant lower limit on |<m>|, typically exceeding 10~ 2 eV, in this case [HHIEE]- 
Using, e.g., the best fit values of \Am\ | and sin 2 one finds: |<m>| ^ 0.02 eV. The 
maximal value of |<m>| is determined by \Am\ | and can reach, as it follows from eqs. 
(0) and (j5J, |<m>| ~ 0.060 eV. The indicated values of |<m>| are within the range of 
sensitivity of the next generation of (/3/?) 0z ,-decay experiments. 

The expression for |<m>|, eq. (}2*2*j) . permits to relate the value of sin 2 a 2 i 1 2 to the 
experimentally measured quantities [2HIEZ1 |<m>|, Am\ and sin 2 29 Q : 

• 2 a 21 ^ ( |<m>| 2 \ 1 

Sm T = V I Ami |cos 4 13 J sin 2 200 ' ( } 

A sufficiently accurate measurement of |<m>| and of |Am^ | and , could allow to get 
information about the value of a^i, provided the neutrino mass spectrum is of the IH type. 

Three Quasi-Degenerate Neutrinos. In this case mo = mi = m 2 = m.3, m\ ^> 
I Am 2 ^ I; m o ~ 0-20 eV. The mass mo effectively coincides with the v e mass m Pe measured 
in the 3 H /3-decay experiments: m = m Pe . Thus, m < 2.3 eV, or if we use a conservative 
cosmological upper limit [12], m o < 0.7 eV. The QD z/-mass spectrum is realized for values 
of mo, which can be measured in the 3 H f3— decay experiment KATRIN [4*T] . 

The effective Majorana mass |<m>| is given by 



< m > = mo 



cos 2 + sin 2 e ia21 ) cos" 0i 3 + e lQ31 sin" 13 (25) 



m cos 2 9 G + sin 2 9 Q e ia21 = m ^1 - sin 2 20 sin 2 (26) 
Similarly to the case of IH spectrum, one has: 

m I cos 200 1 ^ |<m>| ^ m . (27) 

For cos 200 ~ 0.40, favored by the data, one finds a non-trivial lower limit on |<m>|, 
|<m>| ~?t 0.08 eV. For the 90% CL. allowed ranges of values of the parameters one has 
|<m>| ^ 0.06 eV. Using the conservative cosmological upper bound on J^j^j we get 
|<m>| < 0.70 eV. Also in this case one can obtain, in principle, a direct information on one 
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CPV phase from the measurement of |<m>|, m and sin 2 29 Q : 

9 O-zi ( \<m>\ 2 \ \ . . 

sin 2 — = 1 - l - =-!- — T —. 28 

2 V ml j sin 2 26 e v ; 

The specific features of the predictions for |<m>| in the cases of the three types of 
neutrino mass spectrum discussed above are evident in Fig. 2, where the dependence of 
|<m>| on min(rrij) for the LMA solution is shown. If the spectrum is with normal hierarchy, 
|<m>| can lie anywhere between and the presently existing upper limits, eqs. ()12|) and 
(|13|). This conclusion does not change even under the most favorable conditions for the 
determination of |<m>|, namely, even when |Am^ |, Am Q , 6® and #13 are known with 
negligible uncertainty. If the results in [OH] implying |<m>| = (0.1 — 0.9) eV are confirmed, 
this would mean, in particular, that the neutrino mass spectrum is of the QD type. 

5 Implications of Measuring |<m>| 7^ 

If the (/3/3) 0l /-decay of a given nucleus will be observed, it would be possible to determine the 
value of I < m > | from the measurement of the associated half-life of the decay. This would 
require the knowledge of the nuclear matrix element of the process. 

On the NME Uncertainties. At present there exist large uncertainties in the cal- 
culation of the (/3/3)ov-decay nuclear matrix elements (see, e.g., [EH])- This is reflected, in 
particular, in the factor of ~ 3 uncertainty in the upper limit on | < m > \ , which is extracted 
from the experimental lower limits on the (/3/3)oj,-decay half-life of 76 Ge. Recently, encour- 
aging results in what regards the problem of the calculation of the nuclear matrix elements 
have been obtained in [ZH]. The observation of a (/3/3)o y -decay of one nucleus is likely to lead 
to the searches and eventually to observation of the decay of other nuclei. One can expect 
that such a progress, in particular, will help to solve the problem of the sufficiently precise 
calculation of the nuclear matrix elements for the (/3/3) 0l ,-decay [71] . 

Constraining the Lightest Neutrino Mass. As Fig. 2 indicates, a measurement of 
|<m>| ^ 0.01 eV would either i) determine a relatively narrow interval of possible values 
of the lightest v-m&ss min(mj) = m Mm , or ii) would establish an upper limit on m Mm . If 
an upper limit on |<m>| is experimentally obtained below 0.01 eV, this would lead to a 
significant upper limit on m MIN and would imply Am\ > for massive Majorana neutrinos. 



10 



A measurement of |<m>| = {\<m>\) exp ^ 0.02 eV if Am\ = Am^ > 0, and of 



|<m>| = (|< m>\) exp <; \J\Am\ | cos 2 #13 in the case of Am\ = Am\ 2 < 0, for instance, 
would imply that m Mm ^ 0.02 eV and m MIN ^ 0.05 eV, respectively, and thus a i^-mass 
spectrum with partial hierarchy or of QD type [2H] ■ The mass m MIN will be constrained to 
lie in a rather narrow interval [68J (Fig. 2). If the measured value of |<m>|, (\<m>\) exp , 
lies between the min(\<m>\) and max(\<m>\), predicted in the case of IH spectrum, 



\<m> 



± 



Ami I - Ami, cos 2 9 G) ± \l\Am\ | sin 2 9 Q 



cos 2 9, (29) 



m MIN would be limited from above, but min(m MIN )=0 (Fig. 2). If (\<m>\) exp < (\<m>\) max , 
where, e.g., in the case of QD spectrum {\<m>\) max = m M1N = m Pe , this would imply that 
at least one of the two CPV phases is different from zero : «2i 7^ and/or a^i 7^ 10 . 

A measured value of m Pe , (m 9e ) exp > 0.20 eV, satisfying (m 9e ) exp > (m Mm ) max , where 
(m Mm ) max is determined from the upper limit on |<m>| in the case the (/3/3)oi,-decay is not 
observed, might imply that the massive neutrinos are Dirac particles. If (/3/3)ot,-decay has 
been observed and |<m>| measured, the inequality (m Pe ) exp > {m Mm ) max , with (m MIN ) max 
determined from the measured value of |<m>|, would lead to the conclusion that there 
exist contribution(s) to the ((3(3) 0i ,-decay rate other than due to the light Majorana neutrino 
exchange that partially cancels the contribution from the Majorana neutrino exchange. 

Determining the Type of Neutrino Mass Spectrum. The existence of significant 
lower bounds on | < m > | in the cases of IH and QD spectra jSSj , which lie either partially (IH 
spectrum) or completely (QD spectrum) within the range of sensitivity of next generation 
of (/?/3) 0i ,-decay experiments, is one of the most important features of the predictions of 



|<m>|. These minimal values are given, up to small corrections, by \J\Am\ |cos2# Q and 
m cos 29 & . According to the combined analysis of the solar and reactor neutrino data [22J 
i) cos 2^0 = is excluded at ~6<r, ii) the best fit value of cos2# is cos 26*0= 0.44, and iii) 
at 95% C.L. one has for sin 2 13 = (0.02), cos2fl > 0.28 (0.30). The quoted results on 
cos26 l together with the range of possible values of \Am\ \ and m |S1 12H HH EE] lead 
to the conclusion about the existence of significant and robust lower bounds on |<m>| in 
the cases of IH and QD spectrum j^Hl[7S]. At the same time, as Fig. 2 indicates, |<m>| 
does not exceed ~ 0.006 eV for NH spectrum. This implies that max(\<m>\) in the case 
of NH spectrum is considerably smaller than min(|<m>|) for the IH and QD spectrum. 
This opens the possibility of obtaining information about the type of i/-mass spectrum from 
a measurement of |<m>| 7^ 0. In particular, a positive result in the future generation 
of ((3(3) o^-decay experiments showing that |<m>| > 0.02 eV would imply that the NH 
spectrum is strongly disfavored (if not excluded). The uncertainty in the relevant NME and 
prospective experimental errors in the values of the oscillation parameters and in |<m>| 
can weaken, but do not invalidate, these results (see, e.g., ref. 



Constraining the Majorana CPV Phases. The possibility of establishing CP- 
violation due to Majorana CPV phases has been studied in [HE] and in greater detail in 
It was found that it is very challenging 11 : it requires quite accurate measurements 



10 In general, the knowledge of the value of |<m>| alone will not allow to distinguish the case of CP- 
conservation from that of CP- violation. 

^^Pessimistic conclusion about the prospects to establish CP-violation due to Majorana CPV phases from 
a measurement of |<m>| and, e.g., of toq, was reached in . 
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of |<m>| and of mo, and holds only for a limited range of values of the relevant parame- 
ters. For IH and QD spectra, which are of interest, the "just CP-violation" region [28J - an 
experimental point in this region would signal unambiguously CP-violation associated with 
Majorana neutrinos, is larger for smaller values of cos2# Q . More specifically, proving that 
CP-violation associated with Majorana neutrinos takes place requires, in particular, a rela- 
tive experimental error on the measured value of |<m>| smaller than ~ 15%, a "theoretical 
uncertainty" in the value of |<m>| due to an imprecise knowledge of the corresponding 
NME smaller than a factor of 2, a value of tan 2 ^ 0.55, and values of the relevant 
Majorana CPV phases typically within the ranges of ~ (vr/2 — 37r/4) and ~ (5tt/A — 3tt/2). 

6 Conclusions 

Future (/3/3)o u — decay experiments have a remarkable physics potential. They can establish 
the Majorana nature of the neutrinos with definite mass Uj. If the latter are Majorana 
particles, the ov~ decay experiments can provide information on the type of the neutrino 
mass spectrum and on the absolute scale of neutrino masses. They can also provide unique 
information on the Majorana CP-violation phases present in the PMNS neutrino mixing 
matrix. The knowledge of the values of the relevant (/3/3)o„— decay nuclear matrix elements 
with a sufficiently small uncertainty is crucial for obtaining quantitative information on the 
neutrino mass and mixing parameters from a measurement of {(5(3)q v — decay half-life. 

References 

[I] B. T. Cleveland et al, Astrophys. J. 496, 505 (1998); Y. Fukuda et al. [Kamiokande 
Collaboration], Phys. Rev. Lett. 77 (1996) 1683; J. N. Abdurashitov et al. [SAGE Col- 
laboration], J. Exp. Theor. Phys. 95 (2002) 181; T. Kirsten et al. [GALLEX and GNO 
Collaborations], Nucl. Phys. B (Proc. Suppl.) 118 (2003) 33; C. Cattadori, Talk given at 
i/04 International Conference, June 14-19, 2004, Paris, France. 

[2] S. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Lett. B539 (2002) 179. 

[3] Y. Suzuki [Super-Kamiokande Collaboration], these Proceedings. 

[4] Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 87 (2001) 071301; ibid. 89 
(2002) 011301 and 011302; S. N. Ahmed et al, Phys. Rev. Lett. 92 (2004) 181301. 

[5] A. McDonald [SNO Collaboration], these Proceedings. 

[6] Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 81 (1998) 1562. 

[7] K. Eguchi et al. [KamLAND Collaboration], Phys. Rev. Lett. 90 (2003) 021802. 

[8] M. H. Ahn et al. [K2K Collaboration], Phys. Rev. Lett. 90 (2003) 041801. 

[9] B. Pontecorvo, Zh. Eksp. Teor. Fiz. (JETP) 33 (1957) 549 and 34 (1958) 247. 

[10] B. Pontecorvo, Zh. Eksp. Teor. Fiz. 53, 1717 (1967) [Sov. Phys. JETP 26, 984 (1968)]. 

[II] S. M. Bilenky and B. Pontecorvo, Phys. Rep. 41 (1978) 225. 



12 



[12] L. Wolfenstein, Phys. Rev. D17 (1978) 2369; S.P. Mikheev and A.Y. Smirnov, Sov. J. 
Nucl. Phys. 42 (1985) 913. 

[13] S.M. Bilenky and S.T. Petcov, Rev. Mod. Phys. 59 (1987) 671. 

[14] S.T. Petcov, Lecture Notes in Physics 512 (eds. H. Gausterer, C.B. Lang, Springer, 
1998), p. 281; M.C. Gonzalez-Garcia and Y. Nir, Rev. Mod. Phys. 75 (2003) 345. 

[15] A. Bandyopadhyay et al, Phys. Lett. B583 (2004) 134. 

[16] Y. Fukuda et al. [Kamiokande Collaboration], Phys. Rev. Lett. 77 (1996) 1683. 
[17] C. Athanassopoulos et al., Phys. Rev. Lett. 81 (1998) 1774. 

[18] Y. Ashie et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 93 (2004) 101801. 
[19] A. Suzuki [KamLAND Collaboration], these Proceedings. 



[20] E. Aliu et al. [K2K Collaboration], |hep-ex/04"Tl 038, 

[21] V. Gribov and B. Pontecorvo, Phys. Lett. B28 (1969) 493. 

[22] A. Bandyopadhyay et al, Phys. Lett. B608 (2005) 115. 

[23] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28 (1962) 870. 

[24] M. Maltoni et al, |hep-ph/0405172| 



[25] J. Conrad, these Proceedings; M. Sorel, J. Conrad and M. Shaevitz, hep-ph/0305255, 

[26] S.M. Bilenky, J. Hosek and S.T. Petcov, Phys. Lett. B94 (1980) 495. 

[27] J. Schechter and J.W.F. Valle, Phys. Rev. D23 (1980) 2227; M. Doi et al, Phys. Lett. 
B102 (1981) 323; J. Bernabeu and P. Pascual, Nucl. Phys. B228 (1983) 21. 

[28] S.M. Bilenky, S. Pascoli and S.T. Petcov, Phys. Rev. D64 (2001) 053010. 

[29] M. Apollonio et al, Phys. Lett. B466 (1999) 415; F. Boehm et al, Phys. Rev. Lett. 84 
(2000) 3764. 

[30] S.T. Petcov, Phys. Lett. B200 (1988) 373; ibid. B214 (1988) 139 and B406 (1997) 355. 
[31] S.T. Petcov, Phys. Lett. B214 (1988) 259. 

[32] W.C. Haxton, Phys. Rev. Lett. 57 (1986) 1271; S.J. Parke, ibid. 57 (1986) 1275. 

[33] P.I. Krastev and S.T. Petcov, Phys. Lett. B207 (1988) 64; Proc. of the Moriond Work- 
shop on Neutrinos and Exotic Phenomena, Les Arcs, France, January 1988 (ed. J. Tran 
Thanh Van, Editions Frontieres, Gif-sur- Yvette) , p. 173; M. Bruggen, W.C. Haxton and 
Y.-Z. Quian, Phys. Rev. D51 (1995) 4028. 

[34] S.T. Petcov and J. Rich, Phys. Lett. B224 (1989) 401; E. Lisi et al, Phys. Rev. D63 
(2000) 093002. 

13 



[35] P.C. de Holanda, Wei Liao, A. Yu. Smirnov, Nucl. Phys. B702 (2004) 307. 



[36] J.N. Bahcall, M.C. Gonzalez-Garcia and C. Pena-Garay, |hep-ph /0406294 
[37] S.M. Bilenky et al, Phys. Rev. D56 (1996) 4432. 

[38] S. Pascoli and S.T. Petcov, Phys. Lett. B544 (2002) 239; ibid. B580 (2004) 280. 
[39] F. Perrin, Comptes Rendus 197 (1933) 868; E. Fermi, Nuovo Cim. 11 (1934) 1. 
[40] V. Lobashev et al, Nucl. Phys. B (Proc. Suppl.) 91 (2001) 280. 

[41] C. Weinheimer et al, Nucl. Phys. Proc. Suppl. 118 (2003) 279 and these Proceedings. 
[42] M. Tegmark, these Proceedings. 

[43] D.N. Spergel et al. (WMAP Coll.), Astrophys. J. Suppl. 148 (2003) 175. 
[44] M. Lindner, these Proceedings; A. Blondel, these Proceedings. 



[45] M.V. Chizhov, M. Maris and S.T. Petcov, hep-ph/9810501; J. Bernabeu, S. Palomares- 
Ruiz and S.T. Petcov, Nucl. Phys. B 669 (2003) 255; S. Palomares-Ruiz and S.T. Petcov, 
|hep-ph/0406096| 

[46] S.T. Petcov, Phys. Lett. B 434 (1998) 321, (E) ibid. B 444 (1998) 584. 

[47] M.V. Chizhov and S.T. Petcov, Phys. Rev. Lett. 83 (1999) 1096; Phys. Rev. Lett. 85 
(2000) 3979; Phys. Rev. D 63 (2001) 073003. 

[48] V.K. Ermilova et al., Short Notices of the Lebedev Institute 5 (1986) 26; E.Kh. Akhme- 
dov, Yad.Fiz. 47 (1988) 475; P.I. Krastev, A.Yu. Smirnov, Phys. Lett. B226 (1989) 341. 

[49] Q.Y. Liu and A.Yu. Smirnov, Nucl. Phys. B524 (1998) 505; QY. Liu, S.P Mikheyev 
and A.Yu. Smirnov, Phys. Lett. B440 (1998) 319. 



[50] S.T. Petcov, Nucl. Phys. B (Proc. Suppl.) 143 (2005) 159 (hep-ph/0412410). 



[5 1] S. Choubey an d S.T. Petcov, Phys. Lett. B594 (2004) 333; A. Bandyopadhyay et al, 
|hep-ph/0410283; A. Bandyopadhyay et al, Phys. Rev. D67 (2003) 113011. 



[52] P. Langacker et al, Nucl. Phys. B282 (1987) 589. 

[53] V. Barger et al, Phys. Lett. B540 (2002) 247. 

[54] A. de Gouvea, B. Kayser and R. Mohapatra, Phys. Rev. D67 (2003) 053004. 

[55] S. Pascoli, S.T. Petcov and W. Rodejohann, Phys. Lett. B549 (2002) 177. 

[56] S. Pascoli, S.T. Petcov and C.E. Yaguna, Phys. Lett. B564 (2003) 241. 

[57] T. Yanagida, these Proceedings. 

14 



[58] A. Morales and J. Morales, Nucl. Phys. Proc. Suppl. 114 (2003) 141. 
S.R. Elliot and P. Vogel, Annu. Rev. Nucl. Part. Sci. 52 (2002). 



[60] S.T. Petcov, New J. Phys. 6 (2004) 109; S. Pascoli and S.T. Petcov, |hep-ph /0308034 
[61] P. Ramond, these Proceedings; R. Mohapatra, these Proceedings. 



[62] C. Aalseth et al, |hep-ph/041 2300 

[63] L. Wolfenstein, Phys. Lett. B107 (1981) 77; S.M. Bilenky, N.P Nedelcheva and S.T. 
Petcov, Nucl. Phys. B247 (1984) 61; B. Kayser, Phys. Rev. D30 (1984) 1023. 

[64] H.V. Klapdor-Kleingrothaus et al, Nucl. Phys. Proc. Suppl. 100 (2001) 309. 

[65] H.V. Klapdor-Kleingrothaus et al, Phys. Lett. B586 (2004) 198. 

[66] A. Barabash et al, JETP Lett. 80 (2004) 377. 

[67] E. Fiorini, these Proceedings. 

[68] S. Pascoli, S.T. Petcov and L. Wolfenstein, Phys. Lett. B524 (2002) 319; S. Pascoli and 
S.T. Petcov, |hep-ph/0111203| 

[69] S. Pascoli, S.T. Petcov and W. Rodejohann, Phys. Lett. B558 (2003) 141. 

[70] S.M. Bilenky et al, Phys. Lett. B465 (1999) 193; V. Barger and K. Whisnant, 
Phys. Lett. B 456 (1999) 194; F. Vissani, JHEP 06 (1999) 022; M. Czakon et al, 
hep-ph/0003 161} H.V. Klapdor-Kleingrothaus, H. Pas and A.Yu. Smirnov, Phys. Rev. 
D63 (2001) 073005. 

[71] W. Rodejohann, Nucl. Phys. B597 (2001) 110 and |hep-ph/0203214 ; K. Matsuda et al, 
Phys. Rev. D63 (2001) 077301; T. Fukuyama et al, Phys. Rev. D64 (2001) 013001. 

[72] S.T. Petcov and A.Yu. Smirnov, Phys. Lett. B322 (1994) 109. 

[73] V. A. Rodin et al, Phys. Rev. C68 (2003) 044302. 

[74] S.M. Bilenky and S.T. Petcov, hep-ph/ 04052371 

[75] H. Murayama and C. Pena-Garay, Phys. Rev. D69 (2004) 031301. 



15 



